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Abstract The post-Newtonian approximation of general relativistic Liouville's equation is pre- 
sented. Two integrals of it, generalizations of the classical energy and angular momentum, are 
obtained. Polytropic models are constructed as an application. 
Key wordrrelativistic systems:static structures; methods: numerical 

1 Introduction 

Solutions of general relativistic Liouville's equation (grl) in a prescribed space-time have been 
considered by some investigators. Most authors have sought its solutions as functions of the 
constants of motion, generated by Killing vectors of the space-time in question. See for example 
Ehlers (1971), Ray and Zimmerman (1977), Mansouri and Rakei (1988), Ellis, Matraverse and 
Treciokas (1983), Maharaj and Maartens (1985, 1987), Maharaj (1989), and Dehghani and 
Rezania (1996). 

In applications to self gravitating stars and stellar systems, however, one should combine 
Einstein's field equations and grl. The resulting nonlinear equations can be solved in certain 
approximations. Two such methods are available; the post-Newtonian (pn) approximation and 
the weak-field one. In this paper we adopt the first approach to study a self gravitating system 
imbeded in an otherwise flat space-time. In sect. 2, we derive the grl in the pn approximation. 
In sect. 3 we seek static solutions of the post-Newtonian Liouville's equation (pnl). We find two 
integrals of pnl that are the pn generalizations of the energy and angular momentum integrals 
of the classical Liouville's equation. Post-Newtonian polytropes, as simultaneous solutions of 
pnl and Einsteins equation, are discussed and calculated in sect. 4. Section 5 is devoted to 
concluding remarks. 

The main objective of this paper, however, is to set the stage for the second paper in this 
series (Sobouti and Rezania, 1998). There, we study a class of non static oscillatory solutions 
of pnl, different from the conventional p and g modes of the system. They are associated with 
oscillations in the space-time metric, without disturbing the classical equilibrium of the system. 
In this respect they might be akin to the so called gravitational wave modes that some authors 
have advocated to exist in relativistic systems. See, for example, Andersson, Kokkotas and 
Schutz (1995), Baumgarte and Schmidt (1993), Detweiler and Lindblom (1983, 1985), Kojima 
(1988), Kokkotas and Schutz (1986, 1992), Leaver (1985), Leins, Nollert and Soffel (1993), 
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Lindblom, Mendell and Ipser (1997), and Nollert and Schmidt (1992). 



2 Liouville's equation in post-Newtonian approximation 

The one particle distribution function of a gas of collisionless particles with identical masse m, 
in the restricted seven dimensional phase space 

P{m): g^U" = -1 (1) 

satisfies grl: 

£uf = {u*JL - r^U»JL)f(x», TP) = 0, (2) 

where (x^, U' 1 ) is the set of configuration and velocity coordinates in P(m), f(x^, U l ) is a distri- 
bution function, Cu is Liouville's operator in the (x^, U l ) coordinates, and T 1 ^ are Christoffel's 
symbols. Greek indices run from to 3 and Latin indices from 1 to 3. We use the convention 
c = 1 except in numerical calculations of section 4. The four-velocity of the particle and its 
classical velocity are related as 

where f7°(x M ,f*) is to be determined from Eq. (1). In the pn approximation, we need an 
expansion of Ly up to order v 4 , where v is a typical Newtonian speed. To achieve this goal we 
transform (x^,U l ) to (x^,v l ). Liouville's operator transforms as 

where dv 3 jdx^ and dv J / dU % are determined from the inverse of the transformation matrix (see 
appendix A). Thus, 

W U°_.idgap a 



dx» 2Q dx^ 
dvi 1 



v a vP, (5a) 



dU* Q 



v J (goi + 9ikv k ); forz^j, 



(5b) 



-i(^°" 2 + E^(^ + W)); for 

k^i 



where 



Q = U°(g 00 + g ol v l ). (5c) 



Using Eqs.(5) in Eq.(4), one finds 

C v f = U°C v f = 0, 



(6a) 



or 

C v f(x»,v i ) = 0, (6b) 

where 

Cv = V »(A- - ^HvA) - r>^{V ^vi(g 0i + g ik v k )^- 
y dx^ 2Q dx^ dvi' ^ ^ Q dvJ 

d 



~(U°~' + Y, vk (90k +9kiv l )) 7 r-}, (6c) 



We note that the post-Newtonian hydrodynamic equations are obtained from integrations of Eq. 
(6a) over the v-space rather than Eq. (6b) (see appendix B). We expand C v up to the order 
v 4 . For this purpose, we need expansions of Einstein's field equations, the metric tensor, and 
the affine connections up to various orders. Einstein's field equation with harmonic coordinates 
condition, g^T^ v = 0, yields (see Weinberg 1972): 

V 2 2 g 00 = SnG °T 00 , (7a) 
w2 4 _d 2 2 9oo . 2 d 2 2 9oo ,d 2 goo,,d 2 g 00 

-8vrG( 2 T 00 - 2 2 g m °T 00 + 2 T ii ), (7b) 
V 2 3 g i0 = WirG 1 T m , (7c) 
V 2 2 g l3 = —&irG5ij °T 00 . (7d) 

The symbols n g liu and n T liU denote the nth order terms in v in the metric and in the energy- 
momentum tensors, respectively. Solutions of these equations are 

2 goo = -2^, (8a) 
2 g tj = -26 ij( f>, (8b) 

3 9i0 = Vi, ( 8c ) 
4 ff0 o = -20 2 - 2^, (8d) 



where 



0(x,t) = —G / - ± ' ; dV, (9a) 



r x T lU (V t) 

r/'(x,t) = -AG J j x _ x/ | d 3 x', (9b) 

+G 2 T u (x',t)), (9d) 

where bold characters denote the three vectors. Substituting Eqs. (8) and (9) in (6c) gives 

C v = C d + C pn 

d_ j d dcf) d 
dt dx l dx l dv l 
r/A i 2\ d(f) dd> ,■ ,84) dtp 

dru_d_rj 1 dr }l d_ 

where C cl and C pn are the classical Liouville operator and its post-Newtonian correction, re- 
spectively. Equation (6b) for the distribution function /(x M , v l ) becomes 

(£ d + £ pn )f(x»,v i ) = 0. (11) 

The three scalar and vector potentials 4>, ijj and 7/ can now be given in terms of the distribution 
function. The energy-momentum tensor in terms of f(x^, U % ) is 

/TW TJ V 
-^-f^lP^dPU, (12) 

where g = det^g^). For various orders of T^ v one finds 

°T 00 (x x ) = J f(x\v i )d 3 v, (13a) 

2 T V) = j{^v 2 + ^x x ))f{x x ,v i )d 3 v, (13b) 

2 T ij (x x ) = f v i v j f(x x ,v i )d 3 v, (13c) 

^(x*) = J v i f(x x ,v i )d 3 v. (13d) 

Substituting Eqs. (13) in (9) gives 

4>(xL,t) = —G [ f ^ ,t,V p dT', (14a) 
J |x — x'| 

r ] (x,t) = -4Gf Vf ^' ) dr> (14b) 
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3 f v> 2 f(x>,t, 
~2 G J |x-x' 



v' 2 /(x',t,v') dr/ 



^/ /( ;' t,v iy / ' i ;;v^, (Me) 



where cfT = d?xd 3 v. Equations (11) and (14) complete the pn order of Liouville's equation for 
self gravitating systems imbeded in a flat space-time. 

3 Integrals of post-Newtonian Liouville's equation 

In an equilibrium state /(x, v) is time-independent. Macroscopic velocities along with the vector 
potential 77 vanish. Equations (10) and (11) reduce to 

(£ d + £ P*)f( x,v) = - |t A) 

One easily verifies that the following, a generalization of the classical energy integral, is a solution 
ofEq. (15) 

E = \ 2 + 4> + 2(f> 2 + ip + const. (16) 

Furthermore, if </>(x) and ^(x) are spherically symmetric, which actually is the case for an 
isolated system in an asymptotically flat space-time, the following generalization of angular 
momentum are also integrals of Eq. (15) 

k = eij k x j v k exp(-<j)), (17) 

where is the Levi-Cevita symbol. Static distribution functions maybe constructed as func- 
tions of E and even functions of Zj. The reason for restriction to even functions of lf n is to insure 
the vanishing of rf, the condition for validity of Eq. (15). 

4 Polytropes in post-Newtonian approximation 

As in classical polytropes we consider the distribution function for a polytrope of index n as 

F n (E) = ^=(-E) n ^ 2 ; for E < 0, 

= for E > 0, (18) 



where a n is a constant. By Eqs. (13) the corresponding orders of the energy-momentum tensor 
are 



where 



°T n 00 = a n P n (-U) n , (19a) 

2 T°° = a n (3 n <P(-U) n + anln {-U) n+ \ (19b) 

2 C =5 ij 2 T i i = 2a n7n (-Ur + \ (19c) 

V = 0, (19d) 

p n = /V - vT^W^dn = T(3/2)r(n - l/2)/r(n + 1), (20) 
Jo 

7n = /V " r,) n ~ 3/ W /2 dv = r(5/2)r(n - l/2)/r(n + 2), (21) 



and U = 4> + 2(j) + ^ is the gravitational potential in pn order. It will be chosen zero at the 
surface of the stellar configuration. With this choice, the escape velocity v e = ^—2U will mean 
escape to the boundary of the system rather than to infinity. Einstein's equations, Eqs. (7), (8) 
and (9), lead to 

V 2 4> = AttG °T 00 = ATrGa n p n (-U) n , (22) 
VV = 47rG( 2 T 00 + 2 T u ) = 4TrGa n f3 n (f>(-U) n 

+12nGa nln {-U) n+1 . (23) 



Expanding (— U) n as 



(-Ur = (-4>r[l + n(24>+t) + ...], (24) 



ans Substituting it in Eqs. (22) and (23) gives 



V 2 </> = ^Ga n (3 n [{-4>) n - 2n{-<j ) ) n+1 - n(-0) n "V + ••■], (25) 



V 2 ^ = ^Ga n (i n {(^ - \){-4>) n+1 

Pn 

-[3(n + 1)£ - n][2(-4>) n+2 + (-<£)>] +•••}, (26) 

Pn 

These equations can be solved numerically by an iterative scheme. We introduce three dime- 
sionless quantities 

-<f> = X6, (27a) 



-iff = A 2 9, (27b) 
r = < (27c) 

where, in terms of p c , the central density, A = {p c / Oi n [i n ) l ^ n and a~ 2 = 4irGp c /\. Equations 
(25) and (26) in various iteration orders reduce to 

V 2 C O + d n = 0, (28a) 
V^6 + (3^ - 1)^ +1 = 0, (28b) 

Pn 

V^i + 9? = qn{2e n Q +l - 9 n - l ® ), (28c) 

Pn 

g[3(n + 1)£ - n](2C +2 - 0?6 O ), (28d) 

Pn 

where V 2 = ^i^(C 2 -^)- The subscripts and 1 refer to orders of the iteration. The dimension- 
less parameter q is defined as 

_ ^Gp c J _ R s 1 

q ~ # ~ i?2C 1 |^(Ci)|' ( 9) 

where is the Schwarzschild radius, R = a(\ is the radius of system, £i is defined by o (Ci) = 
and c is the light speed. The order of magnitude of q varies from 10 -5 for white dwarfs to 1CT 1 
for neutron stars. For future reference, let us also note that 

-U = X[0 1 + q(@ 1 -2ef)]. (30) 

We use a forth-order Runge-Kutta method to find numerical solutions of the four coupled non- 
linear differential Eqs. (28). At the center we adopt 

HO 

O (O) = 1; #a(0) = ^|o=0; a = 0,1. (31) 

In tables 1 and 2, we summarize the numerical results for the Newtonian and post-Newtonian 
polytropes for different polytropic indices and q values. The pn corrections tend to reduce the 
radius of the polytrope. The higher the polytropic index the smaller this radius. The same is 
true, of course, for higher values of q. 
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5 Concluding remarks 



As we discussed in section 1, some authors is debated to exist a new modes of oscillations in 
relativistic stellar systems. They believed that these modes generated by the perturbations in 
space time metric and no have analogue in Newtonian systems. They used the general rel- 
ativistic hydrodynamics to distinguish them. Although this way is routine but one needs to 
assume some thermodynamic concepts that may be fault in relativistic regime. Hence, to reject 
these conceptual problems, we choosed general relativistic Liouville's equation that is the purely 
dynamical theory. The combination Liouville and Einstein equations enable one to study the 
behavior of relativistic systems without ambiguity. Therefore, in this paper, we used the pn 
approximation to obtain the Einstein-Liouville equation for a relativistic self gravitating stellar 
system. We found two integrals, generalization of the classical energy and angular momentum, 
that are satisfying pnl in equilibrium state. These solutions enable one to construct an equilib- 
rium model for the system in pn approximation. Polytropic models, the most familiar stellar 
models, are constructed in pn approximation. In tables 1 and 2, we compared these models with 
its Newtonian correspondence. The pn corrections tend to reduce the radius of polytrope. The 
higher the polytropic index the smaller this radius. We introduced a parameter q, Eq. (29), to 
enter the effect of central density of system in calculations. Increasing values of q reduce the 
radius of system. In the second paper (Sobouti and Rezania 1998), we study time-dependent 
oscillations of a relativistic system in pn approximation. 
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Appendix A: Derivation of Eqs. (5) 

Consider a general coordinate transformation (X, U) to (Y, V). The corresponding partial deriva- 
tives transform as 





(A.l) 



' dY/dX dV/dX \ I d/dY 
K dY/dU dV/dU ) \ d/dV 

where M is the 7x7 Jacobian matrix of transformation. Setting X = Y = x^, V = v l and 
U = U l for our problem, one finds 

/ 



and its inverse 



One easily finds 



M = 



dx^/dx u dv l /dx v 
dx»/dU j dv l /dU j 



M~ 



dx»/dx v dlP/dx 
dx»/dvi dlP/dvi 



(A.2a) 



(A.2b) 



daf/dx" = 8, 



dx"/dv j = 0, 



dir/dx" = v i du°/dx v = u °y d y v ^, 

' 1 2 dx u 

dU i /dv j = [/%■ + v i dU°/dv j = U 5ij + 17° V^w". 



(A.3a) 
(A.3b) 
(A.3c) 



Inserting the latter in M 1 and inverting the result one arrives at M from which Eqs. (5) can 
be read out. 
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Appendix B: Post-Newtonian hydrodynamics 

Mathematical manipulations in the development of this work has been tasking. To ensure 
that no error has crept in the course of calculations we have tried to infer the post-Newtonian 
hydrodynamical equations from the post-Newtonian Liouville equation derived earlier. From 
Eq. (6a) one has 

CPJ 1 / =U°{C d + C pn )f 

= [(l-0 + iv 2 )£ d + £H/ ; (B.1) 
where C cl and C pn are given by Eq. (10). We integrate C^f over the v-space: 

J CfffjPv = J [(1 - </> + \v 2 )C d + £T]f<Pv. (B.2) 
Using Eqs. (12) and (13), one finds the continuity equation 



T 00 i 2 T 0CK i d , 1 j 3 T 0j\ _ T 00^ _ n 

dt { + > + dxi [ + ' dt ' 



(B.3) 



which is the pn expansion of the continuity equation 

T% = 0, (B.4) 

Next, we multiply CfJ 1 / by v l and integrate over the v-space: 

J tfCfffcPv = j - <j) + ^v 2 )£ d + £ pn ]fd 3 v. (B.5) 

After some calculations one finds 

8 f> 

( lj>0i _|_ 3j>0i\ _|_ ___( 2r £ % 3 _|_ 

dt dxi 

The latter is the correct pn expansion of 

r% = 0; i = 1,2,3. (B.7) 

See Weinberg 1972, QED. 
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Table 1. A comparison of the Newtonian and post-Newtonian polytropes at certain selected 
radii for n=l, 2 and 3 and different values q. 

Table 2. Same as Table 1. n=4 and 4.5. 
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Table 1. 



n 


Polytropic 


Newtonian 


pn polytrope, 6 + q(Q — 29 2 ) 




radius, £ 


polytrope, 9 


q = 10~ 5 


q = 10~ 3 


q = W- 1 




0.0000000 


1.00000 


1.00000 


1.00000 


1.00000 




1.0000000 


0.84145 


0.84145 


0.84143 


0.83936 




2.0000000 


0.45458 


0.45458 


0.45433 


0.42949 


1 


2.9233000 


0.07408 


0.07407 


0.07334 


0.00000 




3.1388500 


0.00087 


0.00086 


0.00000 






3.1415500 


0.00001 


0.00000 








3.1415930 


0.00000 










0.0000000 


1.00000 


1.00000 


1.00000 


1.00000 




1.0000000 


0.84868 


0.84868 


0.84863 


0.84394 




2.0000000 


0.52989 


0.52988 


0.52945 


0.48609 




3.0000000 


0.24188 


0.24187 


0.24031 


0.13289 


2 


3.4737000 


0.13904 


0.13902 


0.13770 


0.00000 




4.3394800 


0.00171 


0.00169 


0.00000 






4.3527000 


0.00002 


0.00000 








4.3529000 


0.00000 










0.0000000 


1.00000 


1.00000 


1.00000 


1.00000 




1.0000000 


0.85480 


0.85480 


0.85473 


0.84773 




2.0000000 


0.58284 


0.58283 


0.58230 


0.52894 




3.0000000 


0.35939 


0.35938 


0.35824 


0.24016 




4.0000000 


0.20927 


0.20925 


0.20764 


0.03226 


3 


4.1939500 


0.18690 


0.18688 


0.18520 


0.00000 




6.8435000 


0.00228 


0.00225 


0.00000 






6.8963000 


0.00002 


0.00000 








6.8967000 


0.00000 
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Table 2. 



n 


Polytropic 


Newtonian 


pn polytrope, 9 + q(@ — 29 2 ) 




radius, ( 


polytrope, 9 


q = 10~ 5 


q = 10~ 3 


q = lO" 1 




0.0000000 


1.00000 


1.00000 


1.00000 


1.00000 




2.0000000 


0.62294 


0.62293 


0.62235 


0.56326 




4.0000000 


0.31804 


0.31802 


0.31645 


0.14194 




5.1541000 


0.22574 


0.22572 


0.22383 


0.00000 




8.0000000 


0.10450 


0.10448 


0.10221 




4 


12.0000000 


0.02972 


0.02970 


0.02716 






14.0000000 


0.00833 


0.00830 


0.00570 






14.6468000 


0.00265 


0.00262 


0.00000 






14.9680000 


0.00003 


0.00000 








14.9713400 


0.00000 










0.0000000 


1.00000 


1.00000 


1.00000 


1.00000 




2.0000000 


0.63965 


0.63964 


0.63905 


0.57857 




4.0000000 


0.36053 


0.36651 


0.35897 


0.18628 




5.7468600 


0.24334 


0.24332 


0.24135 


0.00000 


4.5 


8.0000000 


0.16173 


0.16171 


0.15946 






12.0000000 


0.09015 


0.09013 


0.08766 






16.0000000 


0.05402 


0.05399 


0.05141 






20.0000000 


0.03230 


0.03227 


0.02962 






24.0000000 


0.01782 


0.01779 


0.01510 






28.0000000 


0.00747 


0.00744 


0.00472 






30.2689200 


0.00282 


0.00279 


0.00000 






31.7792300 


0.00004 


0.00000 








31.7878400 


0.00000 
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